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We investigate half-metallicity in [001] stacked (CrAs) n /(GaAs) n heterostructures with n < 3 
by means of a combined many-body and electronic structure calculation. Interface states in the 
presence of strong electronic correlations are discussed for the case n = 1. For n — 2,3 our results 
indicate that the minority spin half-metallic gap is suppressed by local correlations at finite tem- 
peratures, and continuously shrinks upon increasing the heterostructure period. Although around 
room temperature the magnetization of the heterostructure deviates by only 2% from the ideal 
integer value, finite temperature polarization at Ef is reduced by at least 25%. Below the Fermi 
level the minority spin highest valence states are found to localize more on the GaAs layers while 
lowest conduction states have a many-body origin. Our results, therefore, suggest that in these 
heterostructures holes and electrons remain separated among different layers. 
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I. INTRODUCTION 

In recent years, the increasing ability to control the 
growth of semiconductor crystals has made possible the 
fabrication of high-quality artificial heterostructures of 
many different geometries and semiconductor classes. 
Heterostructures interfacing half-metals with semicon- 
ductors are technologically very attractive, since in prin- 
ciple they can be used to attain high polarization spin 
injection from a ferromagnetic electrode into the semi- 
conductor. The high polarization results from the main 
property of half-metallic ferromagnets, namely, the fact 
that they exhibit a metallic density of states for one spin 
channel and a gap at the Fermi-level for the other— 
The major players in the present semiconductor-based 
electronic technology are zinc-blende structures. For this 
reason, half-metals which can adapt to the structure and 
bonding of zinc-blende semiconductors are especially at- 
tractive as they are compatible with existing technology. 

Recently, Zhao and Zunger 3 investigated the relative 
stability of NiAs and zinc-blende (ZB) structures under 
pseudomorphic epitaxial conditions. They found that un- 
der epitaxial growth condition, most of the Cr and Mn 
pnictides and chalcogenides cannot be stabilized below a 
lattice constant of 6.5A However, these conclusions are 
valid for the growth of thick layers, whereas the growth 
of very thin films is dominated by the strain energy at 
the interface. This explains the experimental observation 
that CrAs can be grown in the zinc-blende structure for 
thicknesses up to a few monolayers. 

The first experimental realization of such thin films 
was done by Akinaga et alA. In their work they have 
synthesized zinc-blende CrAs thin films of 3nm thick- 
ness on GaAs substrates. They measured a magnetic 
moment of which is in agreement with theoretical 



prediction, and found an experimental Curie tempera- 
ture above 400 K 4 . From the experimental point of view, 
such small thickness makes this material difficult to use 
in practical devices. Therefore, attention was directed 
towards CrAs/GaAs multilayers^. Previous ab-initio cal- 
culations showed high spin polarization through the en- 
tire region of the multilayer in the case of two monolayer 
CrAs and two monolayer of GaAs stacked alternatively 
(CrAs) 2 /(GaAs) 2^. Initially, the produced multilayers 
of zinc-blende CrAs/GaAs grown on GaAs substrates^ 
indicated that the surface and interface of the multi- 
layer were not completely flat. However, it was found 
that the multilayers grow much thicker than pure zinc- 
blende CrAs. Recently, by optimizing the growth tem- 
perature, the quality was improved significantly, and epi- 
taxial growth of zinc-blende multilayer with flat surface 
and interface was achieved^. The magnetization mea- 
surements for a multilayer structure of (CrAs)2/(GaAs)2 
repeated 100 times^ showed a value of 2/ib per formula 
unit, lower than the theoretical prediction of In ad- 
dition, the temperature dependence of magnetization in- 
dicated a ferromagnetic transition temperature of about 
800K. Moreover, it was confirmed by electronic struc- 
ture calculations that the spin polarization is preserved 
throughout the multilayer, and that it is insensitive to 
substitutional disorder between the Cr and Ga sites^. 

It is the purpose of the present paper to investigate ef- 
fects caused by many-body correlations at finite temper- 
atures in the short period (CrAs)„/ (GaAs)„ heterostruc- 
tures. These ab-initio calculations are performed within 
a combined density functional and many-body dynam- 
ical mean field theory. Our results show that correla- 
tions do not affect much magnetisation in these materials, 
while polarisation is strongly suppressed. In particular, 
the minority spin gap contributes in confining the elec- 
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trons and holes in the heterostructure, an effect similar 
to the one discussed for the (GaAs)„/ (AlAs)„superlatices 
where electrons and holes are spatially separated^. 

The paper is organized as follows: the microscopic de- 
scription of electronic states around the half-metallic gap 
in bulk CrAs is briefly summarized in Section [TTJ The ge- 
ometry of multilayers structures which results in lowering 
the symmetry are discussed in Section IIIII Results are 
presented in Section HVl 



II. BAND-EDGE STATES IN HALF-METALLIC 
BULK CRAS 

Before discussing the nature of the half-metallic gap 
in (CrAs)„/(GaAs)„ heterostructure with or without 
electron-electron interactions, let us briefly summarize 
the principal physical factors leading to gap formation in 
the zinc-blcndc bulk materials. 



states, e g orbitals are practically non-bonding and form 
narrow bands. The bonding / anti-bonding splitting is 
a characteristic of the tetrahedral coordination. The 
Fermi-level is situated in the gap between the bond- 
ing bands and the narrow non-bonding e s -bands. In ad- 
dition, the existence of the band gap is assisted by the 
exchange splitting which keeps minority e g states higher 
in energy. A preliminary analysis of the electronic bands 
in the minority spin channel shows that the nature of the 
gap can change depending on the lattice parameter. For 
a=5.75A, Fig. [5] shows a "direct gap" formed at the X- 
point, while for a=6.06A the gap appears to have an "in- 
direct" nature. In addition, we verified that states at the 
top of the valence band are p-d hybridized states, while 
states at the bottom of the conduction band (electrons) 
are predominantly Cr-d states. As one can see in Fig. 
[2] the entire manifold around the Fermi level is shifted 
to higher energies as the lattice parameter is increased. 
Such changes in the lattice parameter may occur when 
the half-metal is grown on a semiconductor substrate. 




FIG. 1: (Color online) Local-Spin Density Approximation 
(LSDA) - total (black line) and atom resolved (Cr(dashed red) 
As(doted dashed blue)) density of states for bulk CrAs with 
a lattice parameter a=5.75j4 (central panel). Left and right 
panels display the majority and minority spin band structure, 
respectively. States with t2 g and e g symmetry are indicated 
by an arrow at the F point. 



FIG. 2: (Color online) Minority spin bulk- CrAs band struc- 
ture for the different lattice parameters: a — 5. 75 A corre- 
sponding to an enlarged CrAs (long dashed-red), and a = 
6.06A corresponding to the InAs lattice parameter (solid- 
blue). By increasing the lattice parameter the gap changes 
from a direct to an indirect one. Results obtained from a 
LSDA calculation. 



Bulk properties of zinc-blende pnictides and chalco- 
genides were discussed in many papers*^—. In particu- 
lar, first-principles calculation a 4 ' 13 predicted CrAs to be 
half metallic. In this structure every atom has a tetra- 
hedral coordination with the first neighbors being of the 
other atomic species, d— states split into the t2 g and e g 
manifolds. While ^-states hybridize with the p-states of 
the neighboring atom, forming bonding and anti-bonding 



Electronic structure calculations for bulk CrAs for dif- 
ferent lattice constants at finite temperatures and includ- 
ing correlation effects were performed recently by one 
of the authors^. It was shown that correlations induce 
spectral weight in the minority spin gap, while the ma- 
terial remains half-metallic. The spectral weight in the 
minority spin gap is known as so-called non-quasiparticle 
(NQP) states. The occurrence of these states is con- 
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nected to "spin-polaron" processes: the spin-down low- 
energy electron excitations, which are forbidden for half- 
metallic ferromagnets in the one-particle picture, turn 
out to be possible as superpositions of spin-up electron 
excitations and virtual magnon a 15 i 16 . It is important to 
mention that the only situation in which half-metallicity 
is preserved in the presence of finite temperature corre- 
lations is in CrAs and VAs, while in other (semi-)Heusler 
materials the additional spectral weight builds up at 
the Fermi energy, so that strong depolarization takes 
placed—. Depolarisation can originate from other ef- 
fects, such as spin-orbit coupling, which causes a mixing 
of the two spin channels. However, in CrAs this effect was 
found to be less than 1%2£. Consequently, this interac- 
tion is not taken into account for the present calculations. 

There is enormous literature on GaAs band-structure 
calculations, all results showing a direct gap at the T 
point in the band structure. The general discussion 
is focused on the difficulty of density-functional the- 
ory (DFT) mean-field type calculations based on Local- 
Density Approximation (LDA)/ Generalised- Gradient 
Approach (GGA) to reproduce the magnitude of the ex- 
perimental gap which is around leV. It has been as- 
sumed that the band gap problem does not occur in half- 
metals since their dielectric response is that of a metal. 
This assumption was recently contradicted by a GW - 
type of calculation for Lao.7Sro.3Mn03, which predicts a 
half-metallic band gap that is 2 eV larger than the one 
obtained by DFT—. In particular, for the present case 
of semiconducting and half-metal heterostructures this 
might imply that the DFT gaps might underestimate the 
actual experimental values. 

In the present work we investigate the nature of the 
minority spin states involved in the gap formation of 
(CrAs)„/(GaAs)„ heterostructures. In our analysis the 
CrAs bulk states situated at the X-point play a crucial 
role because by Brillouin zone folding associated with 
symmetry lowering present in the (CrAs)„/(GaAs)„ het- 
erostructure, the X point is folded into the T point, 
where the bottom of the conduction band in GaAs is ex- 
pected to be present. Therefore, the relative position of 
the states in X-point with respect to the bottom of the 
GaAs conduction band, would contribute to the band- 
edge states. This effect will be discussed in detail in the 
next sections. Moreover, we investigate up to what ex- 
tent the electronic states forming the (CrAs)„/(GaAs)„ 
band structure are significantly changed by finite tem- 
peratures dynamic correlations. 



III. MULTILAYERS GEOMETRY AND 
METHOD OF CALCULATION 

The (CrAs)„/(GaAs)„ heterostructure, with n < 3 
extends along the [001] direction of the of bulk zinc- 
blende material. The unit cells of the superlattice hav- 
ing the space group symmetry D\ d are simple tetrag- 
onal with c/ag = 2n/v2, where c is the unit-cell di- 



mension along the stacking direction. The basis of the 
tetragonal cell have the constants a = b = a Q / y2, where 
ao represents the FCC lattice parameter. The positions 
of the atoms in the (CrAs) i /(GaAs) i unit supercell are 
Cr:a/2(0,0,0), As:a/2(0, 1, c/a), Ga:a/2(1, 1, c/a) and 
As:a/2(1, 0, 3c/2a) (see Fig. [3]). In the heterostructure 
calculations we included empty spheres with no net nu- 
clear charge in the empty tetrahedral sites in order to ob- 
tain a close-packed structure, as it is generically done for 
zinc-blende semiconductors. For the (CrAs) i /(GaAs) i 
case the empty spheres surrounding the cations and 
anions are located at a/2(l, 0, 3c/2a), a/2(0, 1, c/2a), 
a/2(0, 0, c/a), a/2(l, 1, 0). The positions of the atoms in 
larger unit cells for the other superlattices can be written 
by extrapolating the (CrAs) \j (GaAs) i case. 




FIG. 3: (Color online) (CrAs)\j{GaAs)\ supercell, which 
is repeated twice along the a;-axis for illustrative purposes. 
Atoms situated in different symmetry positions are presented 
by colored spheres Cr (green) - lower and upper planes, As 
(blue) - above and below the Cr layers and Ga (red) - center 
of the unit cell. 

Similar supercells were considered previously in dis- 
cussing the electronic structure of semiconducting 
(GaAs) m /(AlAs)„ superlattices 8 . For the particular case 
when m = n it was demonstrated that within the tetrag- 
onal Brillouin zone the X point (along the k z direction) of 
the FCC zone folds onto the T point. Upon increasing the 
number of layers in the superlattice, the Brillouin zone 
is compressed along k z and the number of zone foldings 
increase. Contrary to the bulk semiconductors where the 
top of the valence band is mainly As-p like states, in the 
case of (GaAs) m /(AlAs)„ heterostructures As-p anionic 
bonds are shared between Ga and Al atoms, leading to 
confinement effects at the valence band maxima. 

For the case of (CrAs)„/(GaAs)„ heterostructures the 
As-p states are shared by Ga and Cr atoms. In par- 
ticular, due to hybridisation the Cr-d states contribute 
significantly to the top of the minority spin valence band. 
Therefore, in addition to the confining effect determined 
by the dimensionality of the heterostructure, electronic 
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correlation are expected to influence the minority spin 
band-edge states around the gap. 

In our calculations we used the lattice parameter 
ao = 5.75A, which is the optimized GaAs lattice con- 
stant obtained from a spin-GGA calculation 22 . This 
value is slightly larger than the one found experimen- 
tally for GaAs (ao = 5.65A) and the one predicted for 
(CrAs) 2 /(GaAs) 2 double monolayers (a = 5.69A)£. We 
checked that our results do not change when the atomic 
sphere radii are chosen to be equal to the average Wigner 
Seitz radius R — 2. 675a. u., in comparison to the case 
were the atomic radii are changed by ±5% depending on 
their type. 

To investigate the effect of electronic correla- 
tions for the above supercells, we performed calcu- 
lations using a recently developed LSDA+ Dynami- 
cal Mean Field (DMFT) scheme 2 ^. Correlation ef- 
fects in the valence Cr-d orbitals are included via 
an on-site electron-electron interaction in the form 
2 Eifm^} ^mm'm"m , "Ci raj c im , (J ,Ci ra '"i,'C im » ff . The in- 
teraction is treated in the framework of DMFT— ~— , 
with a spin-polarized T-matrix Fluctuation Exchange 
(SPTF) type of impurity solver—. Here, Ci ma /c\ mcr de- 
stroys/creates an electron with spin a on orbital m on 
site i. The Coulomb matrix elements t/ mm 'm"m'" can be 
computed for the particular material by taking into ac- 
count the symmetry of the orbitals and the crystal struc- 
ture in terms of effective Slater integrals and Racah or 
Kanamori coefficient a 24 ' 28 . We used the following effec- 
tive Slater parameters: F° = 2eV, F 2 — 5.17eV and 
F 4 = 3.233eV (which results in a Coulomb-interaction of 
U = 2eV and a Hund's rule coupling of J = 0.6eV); these 
are in agreement with previous work o 14 ' 17 ' 29 ' 30 . For the 
exchange correlation functional the LSDA-approximation 
was used, as we found no significant differences with re- 
spect to results using GGA. 

Since static correlations are already included in 
the local spin-density approximation (LSDA), "double 
counted" terms must be subtracted. To obtain this, we 
replace £<x(i?) with T, a (E) — £ CT (0)2i in all equations of 
the LSDA+DMFT procedure^ 4 . Physically, this is re- 
lated to the fact that DMFT only adds dynamical corre- 
lations to the LSDA result. For this reason, it is believed 
that this kind of double-counting subtraction "£(0)" is 
more appropriate for a DMFT treatment of metals than 
the alternative static "Hartree-Fock" (HF) subtraction 3 ^. 

In the calculations we used 287 (6/4) k- vectors for the 
Brillouin-zone integration and used a cutoff of \ m ax—^ 
for the multipole expansion in charge density and the 
potential as well as a cutoff of \ m ax—^ for the wavefunc- 
tions. We checked higher cutoffs for (CrAs)i/(GaAs)i 
and found only negligible differences. 



IV. RESULTS 

A. (CrAs)i /(GaAs)i , heterostructure 

The calculated LSDA-band structure of 
(CrAs)i/(GaAs)i displays an overall metallic be- 
havior. It is of special interest to investigate the 
minority spin channel band structure shown in Fig. 2) 
We plot results in the energy range of —5 to 3eV in order 
to distinguish between the orbitals around the Fermi 
level. In the tetragonal Brillouin zone the symmetry 
points are: T = {w/a)(Q, 0, 0), R = {n/a)(l, 0, a/2), 
A = (7r/a)(l,l,a/c), Z = (tt/o)(0, 0, a/c), 
M = (7r/a)(l,l,0) and X = (7r/a)(l, 0, 0). At 
lower energies (not shown) between —12 and — lOeV one 
can observe a s-like lower valence band of 2eV width. 
This is separated from the upper valence band by an 
interband of 3.7eV. As one can see in Fig. 0]the bands 
in the energy range of —5, —leV are mainly hybridized 
As-p and Cr-d orbitals. At energies around Ep — leV the 
As-p states are hybridized with Cr-d^ and are situated 
at higher energies at the Z-point, with a higher As-p 
weight. In comparison at the top of the valence band in 
the r point a larger weight is visible for Cr-d^ states. 
Here, the As-p contribution is significantly reduced. 
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FIG. 4: (Color online) LSDA minority spin bands for 
(CrAs)i/(GaAs)i along the symmetry-points in the tetrag- 
onal Brillouin zone. Bands are decorated according to dom- 
inant orbitals character. The top of the valence band has 
mainly an As-p character (black), while the bottom of the con- 
duction band has Cr-d^y (red) character. Cr-d^y and d z 2_ 1 
(blue) orbitals cross the Fermi level. 

States around and just above the Fermi level, also show 
interesting features. As we discussed above due to sym- 
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metry reduction, the X-point along k z of the bulk FCC 
zone folds over to the T point of the heterostructure zone. 
Since the lowest conduction band for the minority spin 
spin of bulk CrAs is situated at lower energy than the 
conductoin band edge of bulk GaAs at the T point, one 
would expect that in the (CrAs) 1/ (GaAs) 1 heterostruc- 
ture the lowest states just above the Fermi level are an 
X-like CrAs band. It is clearly seen in Fig. 0] that this 
is not the case, the Fermi level being crossed by Cr-d^y 
as well as d z 2_ 1 orbitals. Note that, for precisely this 
lattice parameter the bulk CrAs and equivalently the 
(CrAs) 1/ (CrAs) 1 is half- metallic. Therefore the pres- 
ence of GaAs determines that Cr-d^y and d z 2_ 1 orbitals 
from the CrAs - which constitutes the "interface-layer" 
- to cross the Fermi level. We have performed selfcon- 
sistent total energy calculations for a fixed volume of the 
unit cell (lattice parameter a = 5.75A) as a function of 
bond length, more precisely as a function of the distance 
between Cr-As and Ga-As atoms along the stacking di- 
rection 6 = zas — z cr which provides information con- 
cerning the stability of the GaAs layers with respect to 
the CrAs ones. In addition this analysis provides also 
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FIG. 5: (Color online) Evolution of LSDA - Cr-A xy minority 
spin density of states of (CrAs)\/(GaAs)\ for different dis- 
tances (S = zas — z cr ) between As and the Cr layer at fixed 
volume of the unit cell. The inset shows a plot of the energy 
(measured with respect to the equilibrium energy E eq ) versus 
S. As one can see, at the optimal 8 = 0.34 the minority spin 
gap is closed and the heterostructure is metallic. 

the evolution of the Cr-d^y states within the minority 
spin gap. As one can see in Fig. [5] within the LSDA 
a larger Cr-As bond length (increasing za s ) and smaller 
GaAs bond favor half-metallicity. The results of the to- 
tal energy calculations show that the equilibrium bond 
length is realized for the distance of <5 ss 0.34 for which 



a metallic (CrAs) 1/ (CrAs) 1 heterostructure is obtained. 
We have checked that LSDA+U calculations with U=2eV 
and J=0.6eV provide the similar value for the equilibrium 
distance. 

Similar calculations were performed including correla- 
tion effects captured by DMFT using two different many- 
body solvers. In fig. [5] we show the comparison between 
the calculations performed using the SPTF and the QMC 
- Hirsch-Fye^ 3 . solvers. The contribution to the total en- 
ergy coming from correlations was computed within the 
LSDA+DMFT-SPTF as an additional Galitskii-Migdal- 
type contribution ^TrE( ? 23 ' 34 , while the interaction en- 
ergy within the LSDA+DMFT-QMC is computed from 
the double occupancy matrix {hi ma hi m i . As one 
can see in Fig. |6] the two numerical results are in good 
agreement. Including correlation effects the equilibrium 




FIG. 6: (Color online) Total energies for different displace- 
ments of GaAs layer with respect the the CrAs layer in the 
n = 1 heterostructure computed within SPTF (blue solid line) 
and QMC (red dashed line) solvers of the DMFT, for U=2eV, 
J=0.6eV. 

distance between the GaAs and the CrAs planes ap- 
proaches the ideal value of 5 = 0.35 (corresponding to 
the FCC structure), thus distortions are not favored in 
this case. Therefore, in what follows we shall consider the 
ideal structure with equidistant GaAs and CrAs planes 
within the supercells. Because of the large supercell 
the following many-body results were obtained using the 
SPTF solver which implements the fully rotational in- 
variant interaction discussed in the previous section. 

In Fig. [7] we compare non-correlated (LSDA) , mean 
field local but static correlations (LSDA+U), with dy- 
namic (LSDA+DMFT) results. All calculations account- 
ing for correlations were performed for the same values of 
the Coulomb parameter U=2eV and J=0.6eV. The mi- 
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to identify the character of non-quasiparticles states. 
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FIG. 7: (Color online) Total and Cr density of states obtained 
by LSDA(green dashed-doted line), LSDA+U (blue dashed 
line) and LDA+DMFT(red solid line) for (CrAs)i/(GaAs)i. 
Results obtained for the ideal S = 0.35. 



nority spin gap opens within the LSDA+U calculation, 
in contrast to the LSDA+DMFT results where spectral 
weight is present in the minority spin gap. As one can see 
in Figj7]in the total density of states (upper panel) the Cr 
states play the dominant role. In the occupied part of the 
majority spin channel, the LSDA states around Ef — leV 
are slightly shifted towards lower energies in LSDA+U, 
while in LSDA+DMFT the same states are shifted to- 
wards the Fermi level, although the spectral weight is 
considerably reduced. At the Fermi level, similar values 
for the density of states are obtained in both mean-field 
calculations, while a slightly larger values are obtained 
including dynamic correlations. For the total density of 
states, Ef ± leV is the energy range were significant dif- 
ferences can be seen. These differences are attributed to 
Cr density of states shown in the lower panel of Fig. [?l 

The LSDA results give a magnetic moment close to an 
integer 2.99/is and by including U at the mean-field level 
an integer magnetic moment of 3/li.b and a gap of 0.94eU 
is obtained. The DMFT density of states is shifted to- 
wards Ef, with states just above the Fermi level. Despite 
their slight spectral weight close to the Fermi level, the 
magnetic moment has a non integer value of 2.88)Ib- We 
attribute this reduction of the magnetic moment to the 
appearance of non-quasiparticle states. To discuss fur- 
ther this effect we plot in Fig. [5] the imaginary/real part 
of minority self-energy for all Cr-d orbitals. One can see 
a strong imaginary part associated with a strong energy 
dependence of the real part predominantly in Ci-d xy and 
Cr-d z 2 orbitals just above the Fermi level which allows 
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FIG. 8: (Color online) Real and imaginary part of the self 
energies for Cr orbitals in the n = 1 heterostructure. Large 
values of 9E for Cr-d x „ (black solid line) and -d z 2_ 1 (dashed 
line) just above the Fermi level are related to the occurrence 
of NQP states. 

The presence of NQP states have been shown previ- 
ously in the bulk CrAs calculations^. Their position 
in bulk is situated at higher energies above the Fermi 
level, in comparison with (CrAs)i/(GaAs)i, so that no 
tails cross the Fermi level and an integer magnetic mo- 
ment of 3(1b is obtained at finite temperatures. From 
the above results we conclude that the presence of GaAs 
layers triggers the presence of "interface" Cr-d xy and Cr- 
d z 2_i states with very small spectral weight within the 
minority spin channel at the Fermi level. In the presence 
of correlations NQP states are formed on these orbitals. 



B. 



(CrAs)„/(GaAs)„ heterostructure, with n=2,3 



For n = 2 the two Cr atoms are equivalent, while the 
As atoms in this heterostructure are shared either by two 
Cr atoms, two Ga atoms or one Cr and one Ga atom. 
In this case the heterostructure can be viewed as two 
(CrAs) i /(GaAs) i interfaces coupled to the GaAs or CrAs 
end of the layers. Increasing the heterostructure further 
to n = 3, makes that Cr atoms become inequivalent, 
having one inner Cr layer and two equivalent outer Cr 
layers. 

For the (CrAs) n / (GaAs)„ heterostructures with n = 
2,3 a self-consistent LSDA calculation yields a half- 
metallic solution with an integer magnetic moment of 
3/Lts per Cr atom. The minority spin gap is 0.94ey for 
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n = 2 and 0.75eV for n = 3. The minority spin gap is 
0.94eV for n = 2 and 0.75eV for n = 3. Within LDA+U, 
the gap is slightly enlarged to l.leV (n — 2), and 0.93eV~ 
(n = 3), while the values for the magnetic moments re- 
main unchanged. 
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FIG. 9: (Color online) Total density of states for 
(CrAs)n/(GaAs) n heterostructures for n = 2 and 3. LSDA 
(green - dashed doted line) and LSDA+U (blue - dashed line) 
calculation yield an half-metallic state. DMFT (red -solid 
line) results suggest the existence of NQP states just above 
the Fermi level with a very small spectral weight at Ef- 

DMFT results for n — 2 show in the majority 
spin channel a redistributions of spectral weight: the 
LSDA/LSDA+U peak at -1/ - l.heV is shifted to lower 
energies, where new structures appear at higher bond- 
ing energies around —5eV. At the Fermi level all three 
methods gives a similar value for the density of states. 
In the minority spin channel NQP states are visible just 
above the Fermi level, while below Ep the band offset is 
similar to the LSDA value. In fig. [9] one can compare 
the overall correlation effects captured at different lev- 
els. Within LSDA+U, the gap is further increased as the 
bandwidth shrinks. On the other hand, dynamic correla- 
tions reduce the gap and enlarge the band width creating 
high-bond energy states and evidence the presence of cor- 
relation induced many-body states in the close vicinity 
of the Fermi level. Similar results are obtained for n = 3. 
The computed values of the magnetic moment are 5.96 
and 8.99/is for n = 2 and n = 3 respectively at 200K. 
At 300K magnetic moments do not change significantly 
Mn=2/3 — 5.97/8.98 and the spin polarization computed 
at Ep is again almost temperature independent and has 
the values of P^ift ~ °- 76 and p dmft ~ °- 71 - 

As a consequence of geometrical relations between the 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 
E-E F (eV) 

FIG. 10: (Color online) Atom and orbital resolved minor- 
ity spin self-energies for (CrAs)2/(GaAs)2 heterostructure for 
U=2eV, J=0.6eV and T=200K. ImE/ReY, are plotted on the 
left /right column. 



heterostructures with n = 1 and n — 2 it is expected that 
also correlation effects would not be very different. In 
he ll 01 we show the layer- and orbital-resolved imaginary 
part of the minority spin self-energy for the heterostruc- 
ture with n = 2. The self-energy plots contain always 
a single curve because of the equivalence of Cr atoms. 
In comparing to Fig. [8] one can see again that the Cr- 
d xy and Cr-d z 2_ 1 states would present similar departures 
from the usual Fermi liquid description, although in this 
case the amplitude of the self-energy anomaly above Ef 
is reduced. This reduction can be attributed to the fact 
that increasing n by adding layers allows for hybridiza- 
tion of CrAs-GaAs layers which softens the self-energy 
anomaly present in the single correlated magnetic Cr 
layer in the n — 1 heterostructure. 

In Fig. [TT] we present the minority spin imaginary 
and real parts of the Crl (inner layer) and Cr2 (inter- 
face layer) selfenergies. To analyze the correlation effects 
taking place for different n, we discuss a comparison of 
the energy dependence of the single particle self-energies 
in Figs. [TOl HU The visible differences are meaningful 
from a physical point of view as a results of the increasing 
number of Cr -layers that can couple with increasing n. 
For n=l the system consists of a correlated single mag- 
netic layer, for which the self-energy apparently deviates 
from the Fermi liquid description. For n — 2 two such 
Cr correlated layers couple and the self-energy anomalies 
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FIG. 11: (Color online) Atom and orbital resolved minor- 
ity spin self-energies for (CrAs)3/(GaAs)3 heterostructure for 
U=2eV, J=0.6eV and T=200K. ImT./Re'E are plotted on the 
left/right column. Crl(interface)/Cr2(inner) atoms are plot- 
ted in dashed-black/red-solid. An almost linear behavior can 
be seen for ReS for inner type atoms. Possible kinks are 
visible for i?e£cvi- 



are softened. Finally, the result for n = 3 does not show 
anomalous features in the energy dependence of the self- 
energy. Upon looking for small energies around the Fermi 
level, one can qualitatively discuss the quasi-particle life- 
time or its inverse the scattering rate 1/rt oc —^s'S(Ep). 
In comparing directly Figs. [T0l [TT] the curvature of Sffi 
is much larger for n = 2 than for n = 3. It is also impor- 
tant to note that for n — 3 the curvature is considerably 
larger for the interface layer, than for the inner layer. 
The above qualitative analysis of the energy dependence 
of the self-energy curves suggest clearly that at the in- 
terface correlation effects are more important leading to 
NQP sates for the n=l and n=2 cases, while for n=3 Cr-d 
states can be described by regular quasiparticle states. 

In order to study the band edges of the minority 
spin channel, we investigate the character (angular mo- 
mentum composition) of the electronic states below and 
above Ep. As discussed above, the bottom of the con- 
duction band consists of NQP states situated just above 
Ep , while the top of the valence band is predominantly 
As-p character. We now discuss the contributions of As- 
p states originating from different As layers to the top of 
minority spin valence band. Fig. Q2] shows the As-p con- 
tributions to the top of the valence band as a function of 
numbers of layers. As the distance between the transition 



metal mono-layer and As mono-layer increases, the As-p 
character continues to increase. At the same time, the 
difference between values from different monolayers de- 
creases. These results suggest that the As layer situated 
closer to the transition metal has its electrons confined 
by the p—d hybridization. At larger distances, hybridiza- 
tion decreases, and more p-character is available to form 
the top of the valence band. At larger temperatures, the 




FIG. 12: (Color online) Relative contribution (in percent) 
of As-p states to the valence band-maxima for each layer 
of the two superlattices (n = 2 and 3). Results are shown 
for (CrAs)n/(GaAs) n heterostructures with U=2eV, J=0.6eV 
and T=200K (red-dashed), 300K (blue-solid). 

As-p relative contribution to the top of the valence band 
slightly decreases within the GaAs-layers and increases 
for As contained within the CrAs layers. This demon- 
strates that the confinement of carriers seen already at 
the LSDA level, i.e. valence-band top almost confined on 
the GaAs region, bottom of conduction band confined on 
CrAs is not significantly changed by finite temperature 
electronic correlations. Therefore the electrons and holes 
remain separated in different regions of the real space. 



V. SUMMARY 

Heterostructures made of layers of Heusler half-metals 
alternated with semiconductors suffer from the presence 
of interface states within the half-metallic gap which 
strongly reduce the possibility of spin-polarized current 
injection. From previous theoretical studies it is known 
that the only case when half-metallicity is retained is the 
case of a NiMnSb(lll)/CdS(lll) interface formed be- 
tween Sb and S atoms^. Alternative options are pro- 
vided by zinc-blende half-metals due to their affinity 
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with the structure of zinc-blende semiconductors. In 
the present paper we have performed combined elec- 
tronic structure and many-body calculations for short 
period (CrAs)„/(GaAs)„ heterostructures, with n < 3 
having a lattice parameter of a = 5.75A. For n = 1 
we investigated the properties of a single-interface layer, 
and demonstrate that Cr-d^ and Cr-d z 2_! states that 
cross the Fermi level acquire a non-quasiparticle charac- 
ter in the presence of correlations. Our results demon- 
strate that heterostructures for n = 2,3 feature mag- 
netizations deviating by less than 2% from the integer 
half-metallic values at temperatures close to room tem- 
perature T = 200, 300K. The computed finite tem- 
perature polarizations deviate significantly from the ex- 
pected 100%. However, we believe that these values 



Pdmft(Ef) « 75%, that are far larger then values for 
other (semi-)Heuslers are still sufficient to make such a 
heterostructure interesting for future spintronic applica- 
tions. In addition, we have investigated the character of 
the minority spin states around the Fermi level. While 
the top of the valence band mainly consists of As-p elec- 
trons localized on GaAs layers, the bottom of the con- 
duction band is mainly determined by many-body in- 
duced NQP states having Cr-d character. This results, 
therefore, suggest that carriers remain separated among 
different layers in these heterostructures. 
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sions. This work was supported by the Austrian science 
fund (FWF project P21289-N16) and by the cooperation 
project "NAWI Graz" (F-NW-515-GASS). 
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